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ABSTRACT 

A new LQG design method is presented which provides prescribed 
imaginary -ax is pole placement for optimal control and estimation systems. 
This procedure contributes another degree of design freedom to flexible 
spacecraft control: Current design methods which interject modal damping 

into the system tend to have little affect on modal frequencies, i.e. they 
predictably shift open-loop plant poles horizontally in the complex plane to 
form the closed-loop controller or estimator pole constellation, but make little 
provision for vertical (imaginary-axis) pole shifts. Imaginary-axis shifts 
which reduce the closed-loop modal frequencies (the bandwidth) are desirable 
since they reduce the sensitivity of the system to noise disturbances. The new 
method drives the closed-loop modal frequencies to predictable (specified) 
levels— frequencies as low as zero rad/sec (real-axis pole placement) can be 
achieved. The design procedure works through rotational and translational 
destabilizations of the plant, and a coupling of two independently - solved 
algebraic Riccati equations through a structured state-weighting matrix. Two 

new concepts, gain transference and Q-equivalency , are introduced and 
employed in the design process. 


1. INTRODUCTION 

Multi-input, multi-output 
systems, such as those encountered 
in flexible spacecraft control, are 
often approached with modern 
optimal control techniques which 
conveniently generate closed-loop 
system gain matrices for simultan- 
eous multi-loop closures. However, 
modern optimal control, as 
presented in most textbooks, is not a 
complete control system design 
methodology. The major problems 
of translating control system 
performance requirements, band- 
width constraints, and compensator 
robustness constraints into the 
performance index have not been 
fully developed [1]. The result is a 
control system design methodology 
that is iterative and empirical. An 
approach to solving these problems 


and de-empiricizing the design 

process is to use structured 
performance index ( S PI ) con- 
straints [2]. SPI constraints may be 
defined as structured performance 
index weighting matrices which 
constrain the weighted variables to 
approach desired predefined 
directions and values in the state 
space as the weighting matrix 
entries approach infinity. This is 
in contrast to generalized con- 
straints for which the weighted 
variables approach zero as the 
weighting matrix entries approach 
infinity. To employ structured 
constraints, and avoid the applica- 
tion of generalized constraints, the 
weighting matrices for the SPI must 
be less than full rank. The poten- 
tial usefulness of the SPI approach 
is apparent: An appropriately- 

structured performance index can 
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drive state variables in predictable 
directions thereby achieving a 
desired performance and bandwidth 
objective. SPFs can provide a non- 
empirical means of constraining 
the controller, estimator, and 
compensator dynamics— the latter is 
critical for closed-loop system 
robustness. 



Figure 1. Design Objective : 
Complex Plane Constraint 


The overall design objective 
for the controller/estimator utiliz- 
ing the SPI approach is illustrated 
in Fig. 1: The gain matrices genera- 

ted through the SPI translate the 
open-loop poles to some prescribed 
closed-loop boundary in the com- 
plex plane. Similarly the compen- 
sation dynamics are constrained to 
reduce closed-loop system sensitivi- 
ty. Note that to achieve this 
objective, two degree-of-freedom 
control is required over each 
controlled mode, i.e. poles require 
movement in two dimensions in the 
complex plane, both horizontally 
(along the real-axis) and vertically 
(along the imaginary-axis). 


In the next section we review 
SPI design methods for prescribed 
real-axis constraint in the optimal 
control, estimation, and 
compensation systems, and 
introduce a coupled Riccati 
equation design technique for 
prescribed imaginary-axis 
constraint. 


2. SPI DESIGN METHODS 

Currently a well-known 
performance index exists for 
prescribed real-axis pole 
translations in the optimal 
controller and estimator systems [1]: 
In the “alpha-shift” technique 
shown in Fig. 2, the standard Linear 
Quadratic Gaussian (LQG) 
performance index is augmented 
with an exponential weighting. 
This exponential weighting 
guarantees that the quadratic terms 
in the performance index decay 
with at least a rate of 2a so that the 
performance index remains finite 
over the infinite interval. The 
result is a guaranteed stability 
margin--all closed-loop poles lie to 
the left of the -2a -line in the 
complex plane. 

The design procedure with 
the alpha-shift technique is 
straight-forward: [+al] is appended 

to the nominal plant dynamics, A. 
This tends to destabilize the plant. 
Optimal control theory is applied 
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and state feedback gains are 
generated for the destabilized plant, 
characterized by [A+al], which are 
guaranteed to stabilize it. When 
these gains are applied to the 
nominal plant, A, the closed-loop 
poles have real parts of -2a. This 
technique provides horizontal 
(real-axis) translation of the plant 
poles from their open to closed-loop 
positions. 

For the compensator, 
predictable real-axis pole 
translations are also possible 
through indirect SPI design 
techniques which structure control 
and observation constraints [3]. 
These constraints tend to normalize 
the control and observation effort 
thereby providing indirect control 
over compensator poles, bandwidth, 
and closed-loop singular values. 

Prescribed imaginary-axis 
pole translations in the optimal 
control and estimation systems are 
the focus of this paper: SPI design 

techniques are presented which 
drive the modal frequencies of the 
closed-loop system to desired levels. 
Conceptually , prescribed 
imaginary-axis pole placement may 
be considered to be composed of a 90 
degree rotation, a vertical 
translation, and a stabilization of 
the open-loop plant poles as shown 
in Fig. 3. Stabilization is achieved 
by generating a stabilization matrix 
for the plant in rotated space and 
applying it to a standard alpha-shift 


design through a SPI. Using the 
stabilization matrix from one 
optimal design process and 
applying it to another couples two 
algebraic Riccati equations (ARE’s) 
together. 

The next section introduces 
two key concepts, gain 
transference and Q-equivalency , 
that are critical to the development 
of the SPI for Riccati equation 
coupling. This is followed by an 
outline of the actual design steps 
required for prescribed imaginary- 
axis pole placement. 


3. DESIGN PROCEDURE: 

PRESCRIBED IMAGINARY- 
AXIS POLE PLACEMENT 

The design procedure for 
prescribed imaginary-axis pole 
placement employs a SPI that 
couples two ARE’s together. Gain 
transference and Q-equivalency are 
important to understanding the 
development of this SPI. 

Gain transference involves 
designing optimal gains for one 
plant and applying them to 
another, indirectly -related, plant. 
As shown in Fig. 4, optimal 
regulator theory is applied to 
system 1, generating optimal gains 

-R’^B^Pj. A closed-loop state feed- 
back system is formed for system 2 
with these gains. (Note that system 



Open-Loop Flint 


1) Rotation 


2) Vertical Translation 3) Stabilization 


Figure 3. Conceptual Development: Prescribed Imaginary-Axis Pole Placement 
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Figure 4. Gain Transference Theory 


2 has identical input vectors as 
system 1.) The optimal gains 
generated for system 1 have been 
transferred to system 2 to form the 
closed-loop system. 

The utility of gain 
transference lies in its harmonic- 
restructuring capability. The 

harmonic structure of the closed- 

1 T 

loop system [A 2 ~BjR Bj P|] can be 

strongly influenced by the 
harmonic structure of Aj. In the 

design procedure for prescribed 
imaginary-axis pole placement, 
optimal gains, Pj, are generated for 

the plant in a rotated space, A^ . 

When these gains are transferred to 
the nominal plant, A 2 , for state 

feedback, the closed-loop system 
takes on the harmonic characteris- 
tics of the plant in a rotated space. 
The state feedback transforms the 
nominal plant to rotational space--a 
key step in achieving prescribed 
imaginary-axis pole placement. 

Q - e q u i v a l e n c y , the other 
concept central to the design 
procedure, involves expanding and 
collecting terms in an ARE to 
indirectly generate a state- 
weighting matrix. An example of 
the concept is shown in Eq. 1 for 
the ARE employed in the alpha- 
shift technique. 


(A+aI) T P+P(A+<xI) - PBR'^P + Q = 0 

A T P + PA - PBR _1 B T P + 2aIP = 0 
Q eq = 2aIP 

= - (A T P+PA-PBR" 1 B T P) (1) 

The alpha terms are 

expanded and collected to form a Q- 

equivalent matrix equal to 2aIP. In 

a SPI, Q is a state-weighting ma- 
eq 

trix that will generate the same 
optimal gains for the nominal 
plant, as those generated through 
the ARE for the alpha-shifted plant. 
This concept is used in the design 
procedure to couple two ARE’s 
together: A Q e matrix for the ARE 

in rotational space is used as the 
state-weighting matrix for an ARE 
in translational (alpha-shifted) 
space . 

An overview of the actual 
design steps that employ the 
concepts of gain transference and 
Q-equivalency are illustrated in Fig. 
5 and described below: 

1 ) Rotational Plant Destabilization. 
A simple matrix transformation 
of the plant rotates poles 
circularly from their open-loop 
positions to the real-axis. This 
removes all harmonic compon- 
ents from the rotated plant 
dynamics. Half of the rotated 
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3) Prescribed Imaginary-Axis 
Pole Placement 


Figure 5. The Three Design Steps for I maginary- Axis Pole Placement 


plant poles in this space are 
unstable. 

2) Rotational Plant Stabilization . A 
stabilization matrix is generated 
for the rotationally-destabilized 
plant through a SPI using 
standard optimal regulator 
design methods. Unstable right- 
hand-plane (RHP) poles at {+CDj, 

+ (D2,..*,+to n } are move d to the left- 

hand-plane (LHP) to positions of 
{-(Op -o> 2 ,...,-co n }, respectively. 

The stabilization matrix does not 
affect the stable LHP poles at {- 
cop -G) 2 ,...,-G) n }. The resulting 

closed-loop system has double 
poles at each modal frequency 
in the LHP. 

3 ) Prescribed Imaginary-Axis Pole 
Placement. The stabilization 
matrix generated for the 
rotationally-destabilized plant is 
used in a SPI to transform an 
alpha-shift design to rotational 
space. The value of a determines 
the closed-loop modal frequen- 
cies, i.e. a prescribes the amount 
of imaginary-axis pole transla- 
tion from the real-axis. 

We now present details of the 
prescribed imaginary-axis pole 
placement design procedure for 
optimal control and estimation 
systems. The three steps outlined 
above are expanded and applied to a 
low-order system to illustrate their 


effects. The complete design 

procedure is then developed and its 
application to flexible spacecraft 
control is illustrated in a numerical 
example. 

DESIGN STEP 1: Rotational 

Plant-Destabilization 

To introduce rotational plant- 
destabilization we compare it 
graphically to the alpha-shift 
technique. As shown in Fig. 6, 
alpha -shifted plant-destabilization 
is accomplished via a horizontal 
translation of the poles into the 
right-half of the complex plane. 
Rotational plant destabilization 
occurs with circular rotations of 
the open-loop poles to the real-axis. 



Figure 6. Plant Destabilization 
Techniques 
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The transformation matrix that 
accomplishes this destabilization of 
the plant is described below for an 
n-mode system. If the plant matrix, 
A, in block-diagonal form is: 



then the rotationally-destabilized 
plant matrix, A r , is defined to be: 

A r = A T 

where 



with eigenvalues X(± 0 )],±( 02 ,...,± 0 ) n ). 
We note that this set of eigenvalues 
is not unique to A r . Other transfor- 
mations of the plant will produce 
block-diagonal or off- diagonal 
matrices with equivalent eigen- 
value sets. These matrices may be 
used in lieu of A r in designing the 
controllers and estimators, but the 
closed-loop system bandwidths tend 
to be larger than those designed 
with A r . A r ’s diagonal structure 
provides the smaller bandwidth 
controller for this design proce- 
dure, and we use this structure in 
the following two-mode example 


illustrating rotational destabiliza- 
tion, and in the flexible spacecraft 
example at end of the paper. 

Given System: to 1 = 10 rad/sec 

(02=15 rad/sec 

One actuator and collocated 
sensor corresponding to torque 
actuation and velocity sensing 
are employed. 



Transformation for A: 



Rotationally-Destabilized Plant: 

10 o o o 

0 10 0 0 

o 0-15 0 

0 O 0 15 

Note that half the eigen- 
values of A r are unstable. Optimal 

regulator design theory may now 
be used to stabilize the rotationally- 
destabilized plant. 


DESIGN STEP 2: Stabilization 

of Rotational Plant 

In this intermediate step, 
optimal regulator theory is applied 
to the plant in rotational space. The 
optimal gains that are generated 
will stabilize the rotationally- 
destabilized plant, but their prac- 
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tical function in this design 
algorithm is to structure a perfor- 
mance index that will rotate and 
stabilize an alpha-shift design for 
the nominal plant. The motivation 
for using the gains in this way 
came from applying them to the 
nominal plant and observing their 
significant effect: They eliminate 

harmonic components from the 
closed-loop system poles. State 
feedback with this set of optimal 
gains rotates the open-loop plant 
poles to the real-axis. This suggests 
that a SPI employing a feedback 
structure with these optimal gains 
as its state-weighting matrix can 
produce rotation and stabilization of 
a prescribed damping design. 

An example of rotational 
stabilization is now presented for 
the low-order system used previous- 
ly. The rotationally-dcstabilized 
plant matrix, A r , developed in step 
1, becomes a parametric design 
matrix in the algebraic Ricatti 
equation (ARE), i.e. the nominal 
plant. A, is replaced with A r in the 
ARE. We note also that the paramet- 
ric design matrix Q is set equal to 
zero in the example. This results in 
double poles in the closed-loop 
systems. Double poles are not 
mandatory. Alternative selections 
of Q may include an identity matrix 
which will split the closed-loop 
system poles, but still maintain 
them in the LHP. Positive scaling of 
the identity matrix will provide as 
much separation of the poles as 
desired. Negative scaling of the 
identity matrix Q adds harmonics to 
the closed-loop system and can be 
used, if desired, to obtain an 
additional increase in the modal 
frequencies in the final design 
step, or to decrease the optimal 
gains. Other structurings of the Q 
matrix are currently being evalua- 
ted for their closed-loop system 
effects. All examples in this paper 
employ a zero matrix Q which 


produces the double-pole structur- 
ing in the intermediate closed-loop 
systems. We now perform the 
rotational stabilization: 


Algebraic Riccati Equation: 

A T P 1 +P 1 A-PjBR' 1 B T P 1 +Q=0 


Parametric Design Matrices: 

A=\ 

R=1 

Q=[01 4x4 

Intermediate Riccati Solution for A r : 


p ,= 


0 0 0 0 

0 5.5556E3 0 -5.0000E3 

0 0 0 0 

0 -5.0000 E3 0 4.6875E3 


Eigenvalues of Intermediate Closed- 
Loop Systems: 

MAh-BR'^Pj) = (-10, -10, -15, -15} 
MA-BR' 1 B T P j ) = (+5, +5, -30, -30} 


Note that Pj is sparse and 
singular. Also, as indicated earlier, 
all harmonic components are com- 
pletely eliminated from the closed- 
loop design model when the optimal 
gains are applied to the nominal 
plant, i.e. all poles have imaginary 
parts equal to zero. The real parts 
are positive or negative values 
which typically have values given 
by one-half or two times the modal 
frequencies. This intermediate 
closed-loop system must now be 
stabilized in a final design step with 
an additional algebraic Ricatti 
cqua-tion which will also add a 
prescribed degree of harmonics to 
the closed-loop system. 


DESIGN STEP 3: Prescribed 

Imaginary-Axis Pole Placement 

In this section we develop 
the SPI that is employed to design 
optimal controllers and estimators 
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with prescribed (closed-loop) modal 
frequencies. The exponentially- 
weighted performance index of the 
alpha-shift technique is modified 
with the optimal gains of the rota- 
tional plant stabilization step. The 
modification of one SPI with the 
optimal gains from another results 
in a coupling of two independently - 
solved ARE's. As developed below, 
the coupling occurs through the 
ARE parametric matrix, Q. 

The optimal gains generated 
in the rotational plant stabilization 
step are used to structure the 
parametric matrix, Q. “Q-equivalen- 
cy”, the expansion and collection of 
terms in an ARE to indirectly 
generate a state-weighting matrix, 
is used to structure Q. A Q C q equa- 
tion, parallel to that shown in Eq. 1, 
is developed for the rotationally- 
destabilized ARE in Eq. 2. (The unity 
subscripts indicate that this is the 
first ARE that is solved in the design 
algorithm.) 

It is Qeql that is used to 
modify the exponentially- weighted 
performance index used in alpha- 
shift designs. The modified perfor- 
mance index and its accompanying 
ARE are shown in Eqs. 3a, 3b, and 4a 


respectively. Qeql transforms the 
alpha-shift design to rotational 
space. After rotation, the alpha 
parameter prescribes the amount of 
imaginary-axis pole translation 
that is desired from the real-axis . 

A Qeq may be developed for 
the modified ARE: Terms in Eq. 4a 

are expanded and collected as 
shown in Eq. 4b. Eq. 4c is formed by 
substitution of Q | and defining 

Qeq2 = 2 aIP 2 . Qeq then is the sum of 
two terms--Q e qi from the rotation- 
ally-destabilized ARE and Qeq2 from 
the alpha-shift design as shown in 
Eq. 5. If Qeq2 » Qeql * i.e. if a is 
large relative to the modal 
frequencies, then the alpha-shift 
term will dominate, and the 
imaginary-parts of the poles will 
asymptotically approach the desired 
alpha value. 

We now demonstrate pre- 
scribed imaginary-axis pole place- 
ment for the low-order system used 
previously. The optimal gains, Pj, 

designed under step 2 are used to 
form Qeql and modify the 
performance index for three alpha- 
shift designs: (ot=0, a=l, and a=2). 


A r Tp i + P,A r - P ] BR 1 B T P 1 = 0 

- PjBR^P, + (A r Pj + PjAj.) = 0 

Q C ql = ( A r P l +P lV = P l BR ' lBTp i < 2 > 

J 2 = J “e 2 “ T [x r Q eq1 x + u 1 Ru]dt (3 a) 

J 2 = J “ e 2 crc [x T P BFT 1 B T FJ x + u T R u) dl (3 b) 

(A+ctI) T P 2 + P 2 (A+al) - P 2 BR' 1 B T P 2 + PjBR^B^Pj = 0 (4a) 

A Tp 2 + P 2 A - P 2 BR 1 B T P 2 + PjBR'^Pj + 2otIP 2 = 0 (4b) 

A T P 2 + P 2 A - P 2 BR l B T P 2 + Qeql + Qeq2 — 0 (4c) 

Qeq = Qeq 1 + Qeq2 = Pj BR ^ B^P j + 2aIP 2 ( 5 ) 


86 



Modified ARE: 


(A+cxI) T P 2 + P 2 (A+aI) - P 2 BR' 1 B T P 2 + PjBR^B^j = 0 


Formed From Optimal Gains P| : 


^eq 1 1 



0 

0 

0 

0 

BR' 1 B T P 1 = 

0 1.1111E5 0 

-1.2500E5 

0 

0 

0 

0 


0 -1.2500E5 0 

1.4062E5 

for Design 1: 

ot=0 




6.2844E5 

-5.1200E4 


-1.0080E6 

3.0400E4 

-5.1200E4 

9.271 1E3 


1.0440E5 

-6.7200E3 

-1.0080E6 

1.0440E5 


1.7601E6 

-6.4800E4 

3.0400E4 

-6.7200E3 


-6.4800E4 

5.3025E3 

for Design 2: 

a= 1 




7.1690E5 

-5.2195E4 


-1.1158E6 

3.0532E4 

-5.2195E4 

1.2205E4 


1.2129E5 

-8.7124E3 

-1.1158E6 

1.2129E5 


1.9794E6 

-7.5562E4 

3.0532E4 

-8.7124E3 


-7.5562E4 

6.6674E3 

for Design 3: 

a=2 




8.2904E5 

-4.8182E4 


-1.2239E6 

2.7139E4 

-4.8182E4 

1.6620E4 


1.401 3E5 

-1.1753E4 

-1.2239E6 

1.4013E5 


2.209 1E6 

-8.7720E4 

2.7139E4 

-1.1753E4 


-8.7720E4 

8.7729E3 


Closed-Loop Eigenvalues ?i(A-BR 


1 B T P 2 ): 


a=0: {-5.0, -5.0, -30.0, -30.0} 

a=l: {-6.1 

a=2: {-7.2 


For a=0, the imaginary-parts 
of the closed-loop poles arc zero 
rad/sec — stable, real-axis pole 
placement is achieved. The inter- 
mediate closed-loop system of step 2, 

characterized by [A-BR" * ], has 

been stabilzed; only the RHP double 
poles at +5 are affected by the new 

- 1 T 

optimal gain matrix, R B P 2 ]. 


± 1.0 i, -31.1 ± 1.0 i} 

± 1.8 i, -32.2 ± 2.0 i} 

For a = l, the imaginary-parts 
of the closed-loop poles have values 
of 1 rad/sec. The real-parts of the 
eigenvalues have been increased 
from their values for the previous 
design as Q e ^2 ^as begun to have an 

effect. 

For a=2, the imaginary-parts 
of the closed-loop poles have values 
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of 2 rad/sec., or approaching 2 
rad/sec. As a increases, Q e q 2 will 

begin to dominate the Q term for 

6 C[ 

the modified ARE--some closed-loop 
modal frequencies may be slightly 
less than the specified modal 
frequencies. 

Step 3 concludes the 
prescribed imaginary-axis pole 
placement procedure for optimal 
controller design. 

The optimal estimator is 
designed via duality theory using 
the same 3 -step procedure. 


We now present the design 
procedure in algorithmic form, and 
illustrate its effects on a higher- 
order example derived from flexible 
spacecraft control. 

4. DESIGN ALGORITHM 

Fig. 7 illustrates the design 
algorithm for prescribed 
imaginary-axis pole placement in 
the optimal controller system. Two 
independently-solved ARE’s are 
employed: the ARE in rotational 

space, and the ARE in translational 
space. The coupling between the 


GIVEN SYSTEM 


x = A x + B u 
y = Cx 



IM (X ( A - B R ’b T P 2 ) ) 


±ai 


1) Transform Plant 
to Rotated Space 


2) Apply Optima! Regu- 
lator Theory to ARE 
in Rotational Space 


3) Form Q-equivalent 
Matrix for ARE in 
Rotational Space 


4) Modify ARE in 
Translational 
Space with Q- 
equi valent Matrix 
From ARE in Rota- 
tional Space and 
Apply Optima! 
Regulator Theory 


5) Form Closed- 
Loop System 


Figure 7. Design Algorithm 
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POSITIVE IMAGINARY-AXIS 


ARE’s occurs with the Q-equivalent 
matrix for the ARE in rotational 
space. This Q-equivalent matrix acts 
as a state-weighting matrix for the 
ARE in translational space. The 
design algorithm for the optimal 
estimator follows a parallel struc- 
ture: Dual variables are substituted 

T 

into ^ and A replaces A in the 
ARE in translational space. 


5. FLEXIBLE SPACECRAFT 
CONTROL EXAMPLE 

The design algorithm is 
applied to a model for the spacecraft 
boom shown in Fig. 8. The model 
contains twelve modes with fre- 
quencies ranging from 0.67 to 11.4 
Hz. Four collocated actuators and 
sensors are positioned at the tip of 
the boom and mid-boom. All modes 
are modeled with zero damping. 

For the example, we design 
five optimal controllers and com- 
pare their pole constellations. The 



Figure 8 . Flexible Spacecraft Boom 

alpha values for the five designs, 
i.e. the prescribed imaginary-axis 
pole placement that is desired, are 
as follows: (cq=0, a2=l» (13=5, 05=10, 
a 5 =15). 

The design results are shown 
in Fig. 9 which plots the 

X(A-BR ^B^P2) for the Five designs. 
(Only the upper-half of the complex 
plane is shown.) Small values of 
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alpha (a 2 = 0 , <*2=1) most closely ap- 
proach the “prescribed” response, 
i.e. the imaginary-parts of all poles 
are approximately equal to the 
prescribed a value. This is due to a 
small or non-existent contribution 

of Q e q2 t0 the Qeq matrix f° r thc 
modified ARE, as explained in the 
previous section. For large values 
of a, the Q e q2 term begins to contri- 
bute to the Q e q matrix, and the 

imaginary-parts of the pole asymp- 
totically approach the desired a 
value. Poles corresponding to 

higher frequency modes have 
imaginary-parts that are closer to 
the a-asymptote. 

6. SUMMARY/FUTURE WORK 

A design procedure has been 

developed for prescribed- 
imaginary axis pole constraints for 
the optimal control and estimation 
systems: The imaginary-parts of 

the closed-loop system poles 
asymptotically approach a 
prescribed value, a. At this stage in 
the development, the maximum 
value that a may assume for a given 
system is constrained, possibly by a 
computational problem with 
solutions for the alpha-shifted ARE, 
Values of a that are large relative to 
the lowest modal frequency in the 
system can produce root migration 
from the desired a-asymptote. 
Small or mid-range frequency 
values of a produce excellent 
results as shown in the example of 
Section 5. Further analysis of the 
computational problem is required. 

The design procedure 
developed empirically as the result 
of numerical experiments in gain 
transference and Q-equi valency 
theory. Future work calls for 
developing an analytical basis for 
the procedure. Additional work 
requires extending the design 
procedure to cover prescribed 


imaginary-axis constraints for the 
optimal compensator system. 
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